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Cells in the animal early visual cortex are sensitive to contour 
orientations and form repeated structures known as orientation 
columns. At the behavioral level, there exist 2 well-known global 
biases in orientation perception (oblique effect and radial bias) in 
both animals and humans. However, their neural bases are still 
under debate. To unveil how these behavioral biases are achieved 
in the early visual cortex, we conducted high-resolution functional 
magnetic resonance imaging experiments with a novel continuous 
and periodic stimulation paradigm. By inserting resting recovery 
periods between successive stimulation periods and introducing a 
pair of orthogonal stimulation conditions that differed by 90° con- 
tinuously, we focused on analyzing a blood oxygenation level- 
dependent response modulated by the change in stimulus orientation 
and reliably extracted orientation preferences of single voxels. We 
found that there are more voxels preferring horizontal and vertical 
orientations, a physiological substrate underlying the oblique effect, 
and that these over-representations of horizontal and vertical orien- 
tations are prevalent in the cortical regions near the horizontal- and 
vertical-meridian representations, a phenomenon related to the radial 
bias. Behaviorally, we also confirmed that there exists perceptual 
superiority for horizontal and vertical orientations around horizontal 
and vertical meridians, respectively. Our results, thus, refined the 
neural mechanisms of these 2 global biases in orientation perception. 

Keywords: fMRI, oblique effect, orientation, primary visual cortex, radial 
bias 

Introduction 

Cells in the early visual cortex in many species of animals are 
sensitive to contour orientations (Hubel and Wiesel 1962, 
1968). Using functional magnetic resonance imaging (fMRI), 
preferences for stimulus orientation have also been inferred 
(Tootell, Hadjikhani, Vanduffel et al. 1998; Boynton and 
Finney 2003; Fang et al. 2005; Liu et al. 2006; Serences et al. 
2009; Weigelt et al. 2012) and decoded (Haynes and Rees 
2005; Kamitani and Tong 2005; Harrison and Tong 2009; 
Swisher et al. 2010) in the human visual cortex. Orientation 
preferences are continuously represented, forming repeated 
structures known as orientation columns (Blasdel and Salama 
1986; Grinvald et al. 1986). High-resolution fMRI has been 
successfully used for visualizing orientation columns in cats 
(Kim et al. 2000; Duong et al. 2001; Fukuda et al. 2006; Moon 
et al. 2007) and humans (Yacoub et al. 2008). 

At the behavioral level, there are 2 well-known global 
biases in orientation perception. One is the oblique effect, re- 
ferring to the observation that both humans and animals have 
a greater behavioral sensitivity to the gratings of cardinal 



(horizontal and vertical) orientations than to the gratings of 
oblique orientations (Appelle 1972). The other is the radial 
bias, describing a biased sensitivity to radial orientations at 
particular angles (Sasaki et al. 2006; Freeman et al. 2011). 
However, it is still unclear how these behavioral superiorities 
for certain orientations are achieved within the early visual 
cortex. Alternative neural mechanisms have been hypoth- 
esized, with one proposing that the more cortical area is 
devoted to achieve this superiority (e.g. Maffei and Campbell 
1970), and the other suggesting that stronger neuronal 
responses could be the source (e.g. Mansfield and Ronner 
1978). In comparison to numerous animal single-unit and 
optical imaging studies, human fMRI studies so far have pro- 
vided only limited evidence that can contribute to this debate, 
as blood oxygenation level-dependent (BOLD) responses to 
different orientations were compared with relatively coarse 
spatial resolutions in the previous studies (e.g. Furmanski and 
Engel 2000; Sasaki et al. 2006; Mannion et al. 2009, 2010; 
Freeman et al. 2011). Thus, it is necessary to address these 
issues at higher spatial resolutions, preferentially at the colum- 
nar (sub-millimeter) resolutions. Furthermore, to deal with the 
increased scan time that accompanies the high-resolution fMRI 
approaches, it is desirable to utilize more efficient stimulation 
paradigms (in comparison with traditional, time-consuming 
block designs) so that orientation preferences at the single- 
voxel level can be revealed within a limited experimental time. 

A continuous stimulation paradigm, in which a grating 
rotates continuously, thus covering all orientations over a 
short period of time, has been used to reveal orientation- 
tuning properties of single voxels in high-resolution fMRI 
studies (Fukuda et al. 2006; Moon et al. 2007; Yacoub et al. 
2008). We have found, however, that the BOLD signal is not 
sensitive enough to capture the modulatory effect caused by 
the change in stimulus orientation when the original para- 
digm is used. The failure is primarily attributed to 2 phenom- 
ena; the first is the gradual decrease in BOLD response, 
possibly caused by neural adaptation mechanisms (Movshon 
and Lennie 1979; Hammond et al. 1985), and the second is 
the phase inconsistency of BOLD response to the repeated 
presentation of particular orientations, noticeably in the latter 
stage of the stimulation period (Bandettini et al. 1997; Mintun 
et al. 2002; Lin et al. 2009). 

To overcome these shortcomings, we devised a stimulation 
paradigm that retains the basic feature of the continuous and 
periodic presentation of the full-range of orientations with 2 
key modifications. First, we inserted a resting recovery period 
following the presentation of the full-range of orientations, 
which allowed a robust BOLD response to be evoked every 
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time when the stimulus appeared, and secondly, we intro- 
duced a pair of orthogonal stimulation conditions, which 
allowed us to extract the maximal BOLD response modulated 
by the change in stimulus orientation using the differential 
analysis method (Grinvald et al. 2000). Using this stimulation 
paradigm, we quantitatively investigated global distribution 
features of single voxels' orientation preferences. We found 
that there are more voxels preferring horizontal and vertical 
orientations at the positions of the horizontal- and vertical- 
meridian representations, thus refining the neural mechan- 
isms of the oblique effect (Appelle 1972) and radial bias 
(Sasaki et al. 2006; Freeman et al. 2011). 

Materials and Methods 

Subjects 

Six subjects (ages 20-35, 3 males) participated in the main fMRI 
experiment. Three of them also participated in 3 additional control 
fMRI experiments and 5 took part in a psychophysical experiment 
measuring discrimination thresholds of line orientations at different 
meridian locations and eccentricities. All subjects were healthy, 
without past history of psychiatric or neurological diseases, had a 
normal or corrected-to-normal vision and gave their written informed 
consent before each experiment. All experiments were approved by 
the RIKEN Functional MRI Safety and Ethics Committee. 



Imaging Hardware 

All fMRI experiments were conducted on a 4 T whole-body MRI 
system equipped with a head gradient system (Agilent Inc., Santa 
Clara, CA, United States of America). High-resolution 3-dimensional 
T,-weighted anatomical magnetic resonance images were acquired 
with a birdcage radio-frequency (RF) coil. A 3-inch single-loop transmit/ 
receive RF surface coil was used to acquire functional and co-registered 
anatomical images in the functional experiment. The surface coil was 
mounted on a support frame attached to the patient table. 

During the experiment, the subject was asked to lie supine on the 
patient table and to rest the back of the head on the surface coil. 
Head motion was restricted using a bite-bar as well as foam rubber 
padding and monitored with 2 pressure sensors placed adjacent to 
the subject's head. The subject's heartbeat was monitored with a 
pulse oximeter, and respiration was monitored with a pressure sensor 
placed on the abdominal region. Both signals were recorded along 
with the timing of RF pulses for later corrections of physiological fluc- 
tuations (Hu et al. 1995). 



Visual Stimulation 

The presentation of visual stimuli in fMRI experiments was controlled 
using the Psychophysics Toolbox extensions (Brainard 1997; Pelli 
1997) in Matlab (The MathWorks, Inc., Natick, MA, United States of 
America) run on a PowerMac G4 computer. Stimuli were delivered to 
the subject's eyes via an optic fiber liquid crystal display goggle 
system (Avotec, Inc., Jensen Beach, FL, United States of America) that 
subtended 25° x 19° of visual angle. The image on the goggles had a 
resolution of 800 x 600 pixels and was presented with a refresh rate 
of 60 Hz. The subject adjusted 2 refractive correction lenses on the 
goggles to achieve a corrected-to-normal vision. Image acquisition 
and stimulus presentation times were synched and a trigger signal 
was sent from the scanner's console to the PowerMac G4 computer 
equipped with a digital I/O board (PCI-DIO-96, National Instruments, 
Austin, TX, United States of America). All signals related to the timing 
of RF pulses, stimulus onset and offset, as well as the state of the 
subject, including physiological data, behavioral performance, head 
position, and eye position, were simultaneously recorded using an 
another PowerMac G4 computer with a data acquisition board 
(PCI-6023E, National Instruments) and a custom-built program 
written in Lab View (National Instruments). 



A circular square-wave grating (18° in diameter, 2 cycles per degree, 
no drifting, 78% contrast) was used in the experiment. In a pilot exper- 
iment, we tested with both stationary and drifting gratings and con- 
firmed that orientation preferences of individual voxels were 
indifferent no matter what type of the grating was used. The stationary 
grating was adopted because the subject reported that stationary 
grating, when rotated in steps (see below), could be perceived more 
smoothly. The foveal region surrounding the fixation cross (up to 1.5° 
from the center of the fixation cross) was filled with homogeneous 
gray color. Throughout the experiment, the subject was instructed to 
fixate on the central fixation cross. To help maintain fixation and the 
level of arousal, the subject performed a detection task on the fixation 
cross. Every 2-5 s, the luminance of the fixation cross was increased 
for 100 ms, and the subject was to report that event with a button 
press. The task itself was only moderately demanding. As a result, all 
subjects performed the task almost perfectly. Nevertheless, the subject's 
eye position was monitored using the eye tracker built into the goggle 
system via the software package iView (SensoMotoric Instruments, 
Boston, MA, United States of America). Offline analysis confirmed that 
the subject maintained a fixation throughout the experiment. 

Imaging Parameters and Experimental Paradigms 

Before the main experiment, the retinotopy of individual subjects was 
mapped using a standard visual-field mapping method (Wandell et al. 
2007). Prior to each functional experiment, Ti-weighted anatomical 
images at the same slice positions as functional images were acquired 
with a 4-segment MPRAGE pulse sequence. In the functional exper- 
iment, 4 contiguous slices (field of view, 16 x 8 cm 2 ; matrix size, 
256 x 128; in-plane resolution, 0.625 x 0.625 mm 2 ; slice thickness, 3 
mm) were placed to cover either the dorsal or the ventral primary 
visual cortex on 1 hemisphere of the subject. This prescription pro- 
cedure, as it has been previously detailed (Cheng et al. 2001), was to 
ensure that a portion of the cortex selected based on morphological 
and retinotopic criteria could be optimally covered with a limited 
number of slices for the sake of a reasonable sampling rate (i.e. 
volume time repetition (TR), see below). 

The stimulus in the main fMRI experiment was continuously and 
periodically presented. The basic idea of the paradigm was to use a 
continuously rotating high-contrast grating to stimulate visually 
responsive voxels. If a voxel is orientation selective, its response 
measured in the BOLD signal should be modulated temporally by the 
rotating grating that changes orientation gradually over the stimu- 
lation time. Further, because the rotation with a full-cycle of orien- 
tations (0-180°) is repeated periodically, the voxel's response should 
also be periodic and a Fourier analysis of the time course at the stimu- 
lation frequency would determine the phase (orientation) at which 
the peak response is evoked, thus allowing us to infer the voxel's 
orientation preference. To exploit the fact that neurons respond maxi- 
mally at their preferred orientations and minimally at orientations 
orthogonal to their preferred orientations (Hubel and Wiesel 1962, 
1968), we adopted the differential analysis method by introducing a 
pair of stimulation conditions whose orientations are orthogonal to 
each other. For the differential BOLD response to this pair of orthog- 
onal conditions, it is expected that the modulatory orientation effect 
(alternating current [AC] components) will be enhanced and other 
response components (direct current [DC] components) that are not 
directly related to orientation changes (e.g. the stimulus onset and 
offset and neural adaptation) will be reduced. 

To focus on the initial stimulation phase, where the problems of 
neural adaptation and loose neurovascular coupling are minimal and 
orientation modulation can be better revealed, we modified the stimu- 
lation paradigm by adding a resting period after each individual 
stimulation period. Specifically, the subject viewed multiple epochs of 
1.5 cycles (36 s) of rotating gratings that were sandwiched between 
24 s blank displays. A full-cycle rotation (0-180°) completed in 24 s 
(this rotation time was empirically chosen based on our preliminary 
experiments where rotation times of 12, 24, and 32 s, were tested), in 
clockwise direction in steps (0.5 s per step; orientation difference 
between successive steps was 3.75°). The 1.5 cycles were chosen to 
compromise between avoiding the sudden strong response at the 
onset of the grating appearance (in case 1 cycle is used) and 
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shortening total scan time (in case 2 or more cycles are used). In indi- 
vidual stimulation periods, the grating started rotating from either the 
horizontal (H) or the vertical (V) orientation, thus creating 2 orien- 
tation conditions that were continuously orthogonal to each 
other. Each orientation condition was repeated 12 times and the 
2 conditions were presented in an H-V-V-H-H-V-V-H-H-V-V-H— 
H-V-V-H-H-V-V-H-H-V-V-H order. Functional images were 
acquired using a high-resolution (0.625 x 0.625 x 3 mm 3 ), 4-segment 
centric-ordered echo-planar imaging (EPI) pulse sequence (4 slices 
per volume; volume TR, 2 s; echo time, 25 ms; average flip angle, 45°) 
to target a retinotopically-defined quadrant of the primary visual 
cortex (VI). In total, 735 volumes of EPI images over 1470 s were col- 
lected in the experiment. 

The same stimulation paradigm was also used in 3 additional 
control experiments to examine the effects of rotation direction, start- 
ing orientations in each pair of orthogonal conditions, and stimulus 
presentation period. In the first control experiment, the rotation direc- 
tion of grating was changed from clockwise to counter-clockwise. In 
the second control experiment, each pair of rotations started from 
either 0° and 90° or 45° and 135°. In the third control experiment, 
the stimulation period of grating was extended from 1.5 to 2.5 cycles. 
The results from these control experiments are summarized in the 
Supplementary Results and Supplementary Figure SI. 



differential time-course (144 volumes). Because the subtraction 
largely removed the DC components that were common between the 
2 orthogonal conditions, this differential time course reflected primar- 
ily the AC component that was modulated by stimulus orientation. 

To determine the significance of orientation modulation contained 
in the differential time course, a fast Fourier analysis was applied and 
the coherence coefficient (cc) was estimated for each voxel. The cc is 
defined as the amplitude of the stimulus frequency (i.e. 12 cycles in 
this case) divided by the sum of squares of all frequency components 
by the following formula (Engel et al. 1997). 

cc= fim . 

^E f \F(f)\ 2 

where / 0 is the stimulus frequency, F is the signal intensity of the Fourier 
transformed time course. Before the Fourier analysis was performed, a 
low-cut de-trending was applied to the differential time course. 

Because the cc reflects the correlation between a sinusoidal func- 
tion and the differential time course (Engel et al. 1997), its significance 
can be estimated, based on the approximation that the cc has an 
F-distribution in the null hypothesis (i.e. zero coherence) with an as- 
sumption of white (temporally uncorrelated) noise, by the following 
formula (Bandettini et al. 1993; Tootell, Hadjikhani, Hall, et al. 1998): 



Data Processing and Analysis 

In all functional scans, longitudinal magnetization was allowed to 
reach steady state before EPI images were collected by the use of 
dummy scans. EPI distortions were minimized using a reference 
volume (without phase encoding) acquired at the beginning of each 
scan (Bruder et al. 1992). The first echo in each segment was a navi- 
gator echo, which was used to correct inter-segment phase and ampli- 
tude variations (Kim et al. 1996). Cardiac and respiratory fluctuations 
were removed from time series images using a retrospective esti- 
mation and correction method, which employed pulsation and respir- 
ation data recorded during image acquisition (Hu et al. 1995). 
In-plane motion correction was also applied (Maas et al. 1997). Both 
corrections were carried out using fe-space data. In addition, a high- 
pass filter was used to suppress baseline signal drifts. No other spatial 
or temporal smoothing was applied. 

For the data obtained in the main experiment, the significance of 
orientation modulation embedded in the raw time course was deter- 
mined on a voxel-by-voxel basis as follows. First, the mean hemody- 
namic response to each stimulus condition was estimated using a 
deconvolution approach (Dale 1999). This procedure assumes a linear 
temporal summation of the fMRI responses (Boynton et al. 1996; Dale 
and Buckner 1997), but does not assume any particular shape for the 
time course of the hemodynamic impulse response function. The 
goodness-of-fit for the deconvolution model was assessed by the 
amount of variance accounted for by the estimated hemodynamic 
responses, r 2 . The significance (P-values) for the r 2 values was com- 
puted using a permutation method described previously (Gardner 
et al. 2005). For each voxel, the P-value gave an estimation of how well 
the voxel responded to the visual stimulation, including both a DC 
component evoked by the high-contrast grating and an AC component 
modulated by the change in stimulus orientation. 

We then created a differential time course for each voxel that 
passed the threshold set for the r 2 value by subtracting neighboring 
pairs of raw BOLD responses to the 2 orthogonal conditions. For each 
of the 12 pairs of H and V conditions as described above, the original 
BOLD response to the H condition was subtracted from that to the V 
condition. The resultant differential time course was thus composed 
of 12 repeats of a 60 s period (30 volumes), with the first 36 s (18 
volumes) exhibiting the difference in response to the 2 orthogonal 
conditions and the remaining 24 s (12 volumes) reflecting the differ- 
ence in the baseline fluctuation. To take the hemodynamic delay into 
account and concentrate on analyzing the differential response modu- 
lated by the change in stimulus orientation, we discarded the first 6 
and last 12 volumes from each period. The retained 12 volumes in the 
period (temporally corresponding to the latter full-cycle stimulation 
with the rotating grating) were then concatenated to form a single 



where n is the number of frequencies. For example, the cc of 0.23 in 
our study (w = 72) corresponds to an f-value of 3. 98, which gives a 
P-value of approximately 0.05. All statistical analyses mentioned above 
were performed using in-house developed functions written in 
Matlab, a Matlab Toolbox for circular statistics (Berens 2009) and the 
software package MrTools developed at the New York University and 
RIKEN Brain Science Institute (http://gru.brain.riken.jp). 

The hemodynamic response delay, which could cause a systemic 
time shift of orientation preference in our study, was corrected by 
comparing the similarity between orientation preferences obtained in 
experiments where the grating was rotated in clockwise and counter- 
clockwise directions, similar to the approach used in the retinotopic 
mapping (Wandell et al. 2007). In brief, the voxels with a significant 
orientation modulation in the clockwise experiment were first se- 
lected using the procedure described above. Then, for the same 
voxels, we obtained differential time courses from the counter- 
clockwise experiment and time-reversed these differential time 
courses. The Fourier analysis was applied to both the differential time 
courses obtained in the clockwise experiment and the reversed differ- 
ential time courses obtained in the counter-clockwise experiment and 
2 sets of phases for all voxels were obtained. We then shifted these 2 
sets of phases in opposite directions (in number of image volumes) 
and calculated the linear correlation of voxel-wise phases between the 
clockwise and counter-clockwise experiments after each shift. The shift 
that gave the largest linear correlation value was defined as the time 
delay in the hemodynamic response. Based on the results from 3 sub- 
jects, the hemodynamic delay in our study was ~ 6 s (3 volumes), 
which was then applied to the data obtained from all 6 subjects. 

Limited by the number of slices that could be prescribed in this 
high-resolution fMRI study (see above), only 1 quadrant of the visual 
field in VI was completely covered for each subject. In total, data 
from 4 dorsal Vis in the right hemisphere, 1 ventral VI in the right 
and 1 ventral VI in the left were collected. In the group analysis, the 
data from different subjects (thus different quadrants) were collapsed 
and represented in a new coordinate system (see Supplementary 
Fig. S2 for explanations). The abscissa of the coordinate represents 
each voxel's retinotopic polar angle, which was converted to within 
0-90° and the ordinate represents the voxel's orientation preference, 
which was converted to take the value either between —90 and 0° or 
between 0 and 90°. In this coordinate system, a voxel at 0° polar 
angle having a 0° orientation preference is a horizontal-meridian 
voxel preferring horizontal orientation, and a voxel at 90° polar angle 
having a —90° or 90° orientation preference is a vertical-meridian 
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voxel preferring vertical orientation; both voxels can be regarded as 
having radial-angle orientation preferences. Similarly, a voxel at 45° 
polar angle having a 45° orientation preference is a voxel having a 
radial-angle orientation preference and a voxel at 45° polar angle 
having a —45° orientation preference is a voxel whose orientation pre- 
ference is complementary to the radial angle. 

Psychophysical Experiment Determining Orientation 
Discrimination Thresholds 

The orientation sensitivity at various locations of the visual field was 
assessed for 5 of the 6 subjects who had participated in the main fMRI 
experiment, using a line orientation discrimination task (Westheimer 
2003). Stimulus locations were distributed along the horizontal (0°/ 
180°), vertical (907270°), and 2 principal oblique (457225° and 1357 
315°) meridians at 2 eccentricities (4° and 8°). The line length was 
0.45° at 4° eccentricity and 0.90° at 8° eccentricity and the width was 
0.04° at both eccentricities. The stimuli were generated using a smooth 
line function in PsychToolbox (DrawLines) on an Apple Mac Pro com- 
puter and presented on a 21" cathode ray tube display (Dell pl013, 
with a resolution of 1800 x 1440 pixels and a refresh rate of 75 Hz) 
through a 2-dimensional Gaussian aperture. The aperture was the same 
as the line length in size and had 2 standard deviations with the full 
transparency centered at the center of the line. The experiment was 
conducted in a darkened room. The luminance of the line and back- 
ground was 83-7 and 0.06 cd/m 2 , respectively. 

At each location, the discrimination thresholds were obtained for 4 
line orientations, namely horizontal, vertical, 45° and 135° orien- 
tations. In a given trial, a reference line for a fixed orientation was 
first presented for 300 ms, followed by a blank screen for 600 ms, and 
then a test line with 1 of 7 orientations (equally spaced intervals, cen- 
tered on the reference orientation) was presented for another 300 ms 
at the same location. After the presentation of the test line, the subject 
was asked to report by pressing 1 of the 2 buttons whether the test 
line appeared to have a more clockwise or counter-clockwise orien- 
tation than the reference line within 2 s. The inter-trial interval was 
3 s. For each reference orientation that was continuously tested, 70 
trials (10 trials for each test orientation) were conducted in a random- 
ized order. In total, 280 trials were collected in a run for 4 reference 
orientations at 1 meridian location. Typically, 4-8 runs were con- 
ducted in an experimental session, and a total of 32 runs (2 runs per 
location) were obtained for all 16 locations at both eccentricities over 
several sessions. For each subject, the test order of reference orien- 
tations and meridian locations at each eccentricity was pseudo- 
randomly assigned. Before the formal experiment was commenced, 
the subject went through sufficient training so that the thresholds 
could be stabilized and proper intervals for the 7 orientations of the 
test line could be empirically determined. 

During the experiment, the subject viewed stimuli binocularly at a 
distance of ~57 cm and was asked to fixate on the central fixation 
cross. The eye position was monitored using an eye tracker (EyeLink 
1000, SR Research, Kanata, Ontario, Canada) and offline analysis con- 
firmed that the subject maintained a fixation during the experiment. 
Data collected across 32 runs were combined and analyzed by the 
method of probits (Finney 1971), using a Generalized Linear Model 
regression function in Matlab. A threshold for each reference orien- 
tation at each meridian location, defined in degrees different from the 
reference orientation, was obtained by taking half the difference 
between the 2 orientations that respectively generated 25% and 75% 
correct responses on the psychometric curve. 



Results 

Revealing Orientation Selectivity of Single Voxels Using 
a Continuous and Periodic Stimulation Paradigm with 
a Differential Analysis Method 

To study orientation-tuning properties of single voxels in the 
human VI, we used a continuous and periodic stimulation 
paradigm (Engel et al. 1994; Kalatsky and Stryker 2003; 



Wandell et al. 2007; Engel 2011) and a differential analysis 
method (Grinvald et al. 2000) (see Materials and Methods). 
The raw time course of a single voxel's BOLD response using 
this stimulation paradigm is shown in Figure 1A (black trace), 
which resembles what is observed typically with a blocked 
design, displaying periodic elevations of BOLD responses in 
between non-stimulated periods. The orientation modulation 
embedded in the raw time course was unveiled as follows. 
First, we used a deconvolution approach (Dale 1999) to 
extract estimated hemodynamic responses to the 2 orthogonal 
gratings, which exhibited typical BOLD response features 
(Fig. IB). After the onset of stimulation, the BOLD response 
rose quickly, yielding a robust sustained (~38-40 s) positive 
response with an overshoot, and after the cessation of stimu- 
lation, the response returned to the baseline gradually. Sitting 
on top of the sustained elevation, there were clear temporal 
modulations seen in both hemodynamic responses that fol- 
lowed the change in stimulus orientation, and importantly, 
the 2 orthogonal conditions resulted in almost opposite 
phases in the modulation (Fig. IB, compare the blue and 
cyan curves in the shaded region). We then obtained a differ- 
ential time course by subtracting the estimated hemodynamic 
response for the horizontal condition from the vertical con- 
dition. This procedure removed larger response components 
that were not directly related to orientation changes, such as 
the sustained response to the presence of a high-contrast 
grating, responses to the stimulus onset and offset, and the 
response decay as a consequence of neural adaptation (DC 
components), and resulted in an enhanced modulatory effect 
that was primarily caused by changes in stimulus orientation 
in the 2 conditions (AC components) (Fig. IB, black trace). 

Quantification of Orientation Selectivity 

To quantify the orientation selectivity revealed by using our 
stimulation paradigm, we assessed the significance of single 
voxels' overall responses to the stimulation and their selectiv- 
ity for stimulus orientation by performing the following 2 ana- 
lyses. We first asked how reliably individual voxels in 
retinotopically-defined regions responded to the visual stimu- 
lation, including the DC elevation of the BOLD response to 
the constant presence of the grating and the AC modulation of 
the BOLD response to the grating rotation. This was deter- 
mined on a voxel-by-voxel basis using a procedure developed 
previously (Gardner et al. 2005). Briefly, we calculated the 
amount of variance in the raw time course (e.g. the black 
trace in Fig. 1A) that was accounted for by the time course 
constructed using the estimated hemodynamic responses (e.g. 
the red trace in Fig. 1A). This value, r 2 , ranging from 0.0 to 
1.0, represents the goodness-of-fit. The smaller the residual 
difference between the raw time course and estimated hemo- 
dynamic responses is, the bigger the r 2 value is. The statistical 
significance of r 2 values was determined in a permutation 
procedure comparing the distribution of another set of r 2 
values calculated using the same data but randomly shuffled 
stimulus times. Because no particular shape for the BOLD 
response was assumed in the deconvolution, both positive 
and negative hemodynamic responses could have survived 
the above procedure. By setting threshold r 2 values CP<0.05) 
for individual subjects, we identified between 5675 and 
12,578 voxels that exhibited significant BOLD responses from 
6 subjects (Table 1, second column). In addition, we also 
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Figure 1. Orientation selectivity revealed by using a continuous and periodic 
stimulation paradigm and a differential analysis method. (4) The raw and estimated 
time courses from a representative V1 voxel. The BOLD response was measured, 
while the subject viewed a grating that started to rotate either in horizontal (indicated 
by cyan arrows) or in vertical (blue arrows) orientation. Each rotation lasted 1.5 
cycles, and critically, was followed by a 24 s blank period. A full-cycle (0-180°) 
rotation took 24 s and there were 12 repeats for each rotation condition in a 
functional scan. The initial orientation and rotation direction of the grating are 
schematically indicated above the time courses and the arrows point at the start 
times of individual rotations. (6) Estimated hemodynamic responses of the voxel to 
the grating that started to rotate in the horizontal (cyan) and vertical (blue) 
orientations, respectively. The black curve shows the difference between the 2 
estimated hemodynamic responses (blue minus cyan curve). The horizontal bar 
indicates the stimulation period and the shaded area indicates the final full-cycle 
stimulation between 12 and 36 s. (C) Time course showing the difference between 
BOLD responses to the grating that started rotation from the vertical orientation (e.g. 
V, and V 2 ) and that to the grating that started rotation from the horizontal orientation 
(e.g. H, and H 2 ) for the first 2 pairs of gratings. Shaded areas indicate differential 
responses in the final full-cycle of rotation between 1 2 and 36 s. (D) The differential 
time course constructed by concatenating 12 truncated differential responses 
obtained in the final full-cycle. The red curve is the sinusoidal fit to the differential 
time course. 



Table 1 

Number of voxels after each selection criterion 



Subject 


Number of voxels (%) 


Number of voxels (%) 


Number of 




after selection with r 2 


after selection with 


voxels (%) 




value, P < 0.05 


estimated hemodynamic 


after selection 






response <0.5% 


with cc value, P 








< 0.05 


SI 


10,864 (100) 


8950 (82.4) 


3955 (36.4) 
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S3 


8731 (100) 


7675 (87.9) 


3612 (41.4) 


S4 


12,578 (100) 


9069 (72.1) 


3794 (30.2) 


S5 


6851 (100) 


3967 (57.9) 


2946 (43.0) 


S6 


5675 (100) 


4569 (80.5) 


1727 (30.4) 


Mean (%) 


100 


76.7 


35.3 
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0.15 


B 0.16 
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Coherence coefficient (cc) Coherence coefficient (cc) 

Figure 2. A voxel's orientation selectivity was assessed by its cc that reflects 
sinusoidal modulation at the stimulation frequency. (4) Distribution of ccs of all voxels 
that responded significantly to the visual stimulation for a representative subject. Left 
arrow indicates the mean of the distribution and right arrow indicates the value of cc 
at P < 0.05. (£?) Averaged distributions of ccs obtained from 6 subjects. 



Table 2 

Number of voxels in the studied quadrant of VI 



Subject 


Quadrant of V1 


Number of voxels (%) 
inside the quadrant 


Number of voxels (%) not 
associated with large surface 
veins (estimated hemodynamic 
response <3%) 


SI 


Right dorsal VI 


661 (100) 


608 (92.0) 


S2 


Right dorsal VI 


402 (100) 


352 (87.6) 


S3 


Right dorsal VI 


499 (100) 


463 (92.8) 


S4 


Right dorsal VI 


616 (100) 


545 (88.5) 


S5 


Right ventral V1 


489 (100) 


489 (100) 


S6 


Left ventral V1 


308 (100) 


294 (95.5) 


Mean (%) 




100 


92.7 



introduced another criterion (estimated hemodynamic 
response >0.5%) typically used in our laboratory for studies 
employing similar visual stimuli to remove the voxels that 
were weak in response and tended to be located outside the 
gray matter (likely due to partial voluming). In total, between 
3967 and 9069 voxels from 6 subjects survived and were re- 
garded as visually responsive voxels in the present study 
(Table 1, third column). 

We then determined the statistical significance of orien- 
tation modulation for the voxels that passed the above 
thresholds. For this purpose, we first constructed a differential 
time course from its original time course (the black trace in 
Fig. 1A) for each voxel. This differential time course was com- 
posed of 12 pair- wise subtractions between neighboring 
BOLD responses to the grating that started rotation from the 
vertical orientation (e.g. Fig. 1A, Vj and V 2 ) and to the grating 
that started rotation from the horizontal orientation (e.g. 



Fig. 1A, Hi and H 2 ). For each of 12 resultant subtracted BOLD 
responses (e.g. Fig. 1C, V x — Hj and V 2 -H 2 ), the original 60 s 
time course was truncated to retain only the 24 s period corre- 
sponding to the final full-cycle (0-180°) orientation change 
(Fig. 1C, shaded ranges). /All 12 truncated responses from the 
voxel were eventually appended to form a single differential 
time course, which reflected primarily the AC component that 
was repeatedly modulated by stimulus orientation (Fig. ID, 
black trace). We applied the Fourier analysis to the differential 
time courses of these selected voxels and calculated their ccs. 
For a representative subject, these values for the voxels that 
responded significantly to overall visual stimulation, ranged 
from 0.01 to 0.75, with a mean of 0.15 (Fig. 2A; across all the 
6 subjects, 0.16 ±0.03, mean ± standard deviation [SD], 
Fig. 2E). Because a cc reflects the correlation between a sinu- 
soidal function of a particular frequency and a voxel's differ- 
ential time course (Engel et al. 1997), we estimated its 
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significance based on the approximation that the cc follows 
an F-distribution in the null assumption, that is, zero coher- 
ence (Bandettini et al. 1993; Tootell, Hadjikhani, Hall, et al. 
1998). We set the threshold for the cc at 0.23 (P<0.05). On 
average, slightly more than a third of the voxels that exhibited 
significant BOLD responses passed this threshold (35-3%) and 
were regarded as orientation selective (see Table 1, fourth 
column, for numbers of selected voxels from individual sub- 
jects). It should be emphasized that the orientation preference 
determined using this approach did not depend on the 
settings of particular parameters (see Supplementary Results 
and Supplementary Fig. SI for the results from control 
experiments). 



Biased Orientation Representations in Human VI 

Using the phase information derived from the Fourier analysis 
to estimate orientation preferences of activated voxels, we 
quantitatively investigated global distribution features of indi- 
vidual voxels' orientation preferences. Because the number of 
slices was limited in this high-resolution fMRI study, the 
results described below were from only 1 quadrant of VI of 
each subject and were generalized across 6 subjects (see 
Table 2 for number of voxels and studied quadrants from indi- 
vidual subjects). 

We first found that there are more voxels preferring hori- 
zontal and vertical (cardinal) orientations. Figure 5A illustrates 
these biases for a representative subject, from whom the data 
were acquired from the right dorsal VI, representing the left 
lower visual field in the parafoveal region (2-9° eccentricity). 
Across all 6 subjects, this tendency was significant (1-way 
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Figure 3. Horizontal and vertical orientations are over-represented in human V1. (4) 
The number of orientation-selective voxels plotted against the preferred orientation 
for a representative subject. (fl) The averaged percent BOLD responses of voxels 
plotted against their preferred orientations. (C) Proportion of orientation-selective 
voxels preferring cardinal (around 0° and 90°) and oblique (around 45° and 135°) 
orientations across 6 subjects. Each bar represents pooled orientation preferences of 
single voxels in a 30° range, centered at 0°, 45°, 90°, or 135°. Voxels with an 
estimated preferred orientation at the borders between cardinal and oblique 
orientations (e.g. 15°-30°) were excluded in the analysis to avoid the possible 
overlap between the populations. (D) The averaged percent BOLD responses of 
voxels tuned to cardinal and oblique orientations, averaged across all 6 subjects. Error 
bars in C-D indicate 1 SD. *P < 0.05. 



analysis of variance (ANOVA), P< 0.001; Fig. 30- Multiple- 
comparison tests also revealed that the voxels preferring car- 
dinal orientations were significantly more (P<0.05) than 
those preferring oblique orientations, whereas there was no 
significant difference between voxels preferring the 2 cardinal 
orientations or between those preferring the 2 oblique orien- 
tations. It should be emphasized that the biased population 
orientation selectivity observed here did not depend on the 
starting orientations in the pair of orthogonal conditions (see 
Supplementary Figs SI C,D). The biased selectivity did not 
depend on which stimulation period was chosen, either, as 
we obtained almost identical results (data not shown) when a 
sliding full-cycle time window was used, which was shifted 
gradually within the 1.5 cycles, after hemodynamic delay was 
adjusted. It has been well documented that both humans and 
animals have a greater behavioral sensitivity to the gratings of 
cardinal orientations than to the gratings of other oblique 
orientations, a phenomenon known as the oblique effect (see 
review in Appelle 1972) and alternative neural mechanisms 
underlying the oblique effect have been hypothesized, with 
one proposing that there are more neurons preferring cardi- 
nal orientations in the visual cortex (Maffei and Campbell 
1970) and the other suggesting that the oblique effect may be 
related to greater responses of neurons tuned to cardinal 
orientations in the human visual cortex (Mansfield and 
Ronner 1978; Furmanski and Engel 2000). In our study, there 
was no evidence for the second mechanism. Although the 
averaged strength of BOLD responses (defined by the magni- 
tude value from the Fourier analysis) of the voxels tuned to 
cardinal orientations tended to be larger than that of the 
voxels tuned to oblique orientations in some subjects 
(Fig. 3B), this trend was not significant across 6 subjects 
(1-way ANOVA, P= 0.41, Fig. 30). 

We then examined the relationship between individual 
voxels' orientation preferences and their retinotopic polar 
angles, which is related to another global orientation bias (the 
radial bias) found in the early visual cortex, including VI (Le- 
venthal 1983; Leventhal et al. 1984; Schall et al. 1986; Sasaki 
et al. 2006; Freeman et al. 2011). Radial orientations are like 
wheel-spokes protruding from the central fixation point. The 
voxels showing the radial bias should have the orientation 
preferences that are consistent with their retinotopic polar 
angles, for example, the voxels located around 45° polar 
angle have dominant preferences for orientations near 45° 
and the voxels located around the horizontal/vertical meri- 
dian representations have preferences near the horizontal/ver- 
tical orientations. This relationship for a representative subject 
is depicted in Figure 4A. Examining Figure 4A, as well as 
Figure 4D, where the data from all 6 subjects are pooled, we 
found that there appeared to be over-representations for 
horizontal orientations (around 0° orientation) near the 
horizontal-meridian representation (around 0° polar angle) 
and vertical orientations (around 90° and —90° orientations) 
near the vertical-meridian representation (90° polar angle). To 
confirm this impression, we analyzed the relationship 
between a voxel's preferred orientation at a given polar angle 
and the radial orientation expected for the polar angle, and 
classified this relationship into 4 orientation categories, that 
is, the preferred orientation either matched the radial orien- 
tation (Radial) or differed from the radial orientation by +45° 
(Radial +45), -45° (Radial -45) or 90° (Tangential). By group- 
ing the data into 3 ranges of polar angles (Fig. 4B for the 
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Figure 4. The radial bias manifested by the relationship between individual voxel orientation preferences and their retinotopic polar angles. \A) Preferred orientations of the 
voxels from the same subject shown in Figure 3A plotted against their retinotopic polar angles, (fl and C) Proportion and percent BOLD responses of the voxels classified into 1 
of 4 categories of relationship between a voxel's preferred orientation at a given polar angle and the radial orientation expected for the polar angle, grouped in 3 ranges of polar 
angles. (D, F, and F) are the same as Ifi, B, and C), respectively, but are from the data pooled from all 6 subjects. Error bars in F and F indicate 1 SD. There were significant 
radial orientation biases for the voxels near horizontal- and vertical-meridians (F, see text for statistical details), but these biases were not significant in terms of percent BOLD 
responses (F). 



representative subject and Fig. 4E for all 6 subjects), a 2-way 
ANOVA revealed a significant main effect (P< 0.0001) of 
orientation category and a significant interaction (P< 0.0001) 
between orientation category and the range of polar angles. 
Subsequent 1-way ANOVA and multiple-comparison tests for 
each of the 3 ranges of polar angles indicated that there were 
over-represented voxels preferring horizontal and vertical 
orientations at horizontal- and vertical-meridian represen- 
tations (P< 0.0005 and P< 0.005, respectively). Away from the 
horizontal and vertical meridians, in the range of polar angles 
between 22.5° and 67.5°, only vertical orientations (Radial 
+45°, Fig. 4E) were found to be over-represented (P< 0.0001). 
There was a slight tendency for the radial bias (Radial > Tan- 
gential), but this difference was not significant CP = 0.95, 
Fig. 4E). These observations, thus, suggest that the radial bias 
indeed exists in VI, but the bias is spatially heterogeneous: It 
is very strong around horizontal and vertical meridians, but 
negligible around oblique meridians. These biases, however, 
were not observed in the strength of BOLD responses (2-way 
ANOVA, ,P=0.82, P=0.69, P=0.24 for the orientation cat- 
egory, range of polar angles and their interaction, respect- 
ively; Figs 4C,F). 

Due to the limitation in brain coverage in this high- 
resolution fMRI experiment, we were not able to systemati- 
cally investigate the dependency of the oblique effect and 
radial bias on spatial locations. It may be argued that these 
biases cannot be generalized across the entire visual field 
because the data were averaged across different visual quad- 
rants; it is likely that some quadrant-dependent discrepancies 
exist (Van Essen et al. 1984). To examine the possible depen- 
dency on the eccentricity (Vandenbussche et al. 1986; Davey 
and Zanker 1998; Westheimer 2003), we have split the data 
into 2 halves with roughly the same number of voxels, and 
performed similar analyses on the 2 datasets. The results 
suggested that the oblique effect and radial bias were not 



different across different eccentricities (see Supplementary 
Figs S3 and S4). 

To exclude the possibility that biased orientation represen- 
tations observed in our study were due to the presence of 
spatially irregularly distributed macrostructures (Boynton 
2005; Swisher et al. 2010), such as large surface veins, which 
may carry orientation-selective information (Gardner 2010; 
Kamitani and Sawahata 2010; Kriegeskorte et al. 2010; Op de 
Beeck 2010; Shmuel et al. 2010), we removed the voxels 
associated possibly with large surface veins by settings an 
upper threshold of the estimated hemodynamic response of 
3%, which could reasonably well separate the voxels located 
in the brain tissue and those in the veins in our studies 
(Cheng et al. 2001; Sun et al. 2007). Across 6 subjects, we 
found that all the features described above retained for the 
voxels that were devoid of large surface veins (92.7% of total 
VI voxels; see Table 2, third column): There were over- 
representations for cardinal orientations (1-way ANOVA, 
P< 0.005) and the voxels preferring cardinal orientations 
were more than those preferring oblique orientations 
(multiple-comparison tests, P<Q.Q5) (Fig. 5A), cardinal orien- 
tations were predominantly represented around respective 
meridian representations (Fig. 5C, see legends for details of 
statistical tests), and these biases were not reflected in the 
BOLD response strength (Fig. 5B,D). 

Biased Orientation Sensitivity Observed at the 
Behavioral Level 

In an attempt to directly relate our findings regarding biased 
orientation representations to possible perceptual superiority 
for these orientations, we conducted a psychophysical exper- 
iment measuring the orientation sensitivity at various meri- 
dian locations (horizontal, vertical, and 2 principal oblique 
meridians) and eccentricities (4° and 8°, which fall within the 
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size of the grating used in the fMRI experiment) (Fig. 6A). 
The 5 observers in this experiment were the subjects who had 
participated in the fMRI experiment. In contrast to the fMRI 
measurement, where individual results were obtained from 
single quadrants, the discrimination thresholds for line orien- 
tations were measured throughout the visual quadrants. Con- 
sistent with the fMRI results, we found a clear oblique effect 
behaviorally in our subjects. Across meridian locations and ec- 
centricities, the discrimination thresholds for horizontal and 
vertical orientations were smaller than those for oblique 
orientations. To compare with the fMRI results directly, we 
analyzed the behavioral data in the same way as we did for 
the fMRI data. Similar to grouping fMRI data into 3 ranges of 
polar angles (Figs 4B,E, and 50, we collapsed behavioral 
data from all measured locations into 3 meridian locations, 
namely the horizontal meridian, vertical meridian, and other 
principal oblique meridians and converted the relationship 



B 



C 7 



I 




0± 15 45 ±15 90 ±151 35 ±15 

Preferred orientation (") 



■ Radial □ Tangential 
□ Radial+45 □ Radial-45 , 



0±15 45±1590±15 135 ±15 

Preferred orientation (°) 




D 



.2 1 1.5 

Q 
_l 
O 




0-22.5 22.5-67.5 67.5- 

Polar angle (°) 



0 - 22.5 22.5 - 67.5 67.5 - 90 
Polar angle (°) 



Figure 5. Preference for cardinal orientations and the radial bias are preserved after 
the removal of the voxels associated with large surface veins. The conventions in (4) 
and (S) are the same as in Figure 3C and D, respectively, and those in (C and D) are 
the same as in Figure 4f and F, respectively. In (C), a 2-way analysis of variance 
(ANOVA) revealed a main effect of orientation category and a significant interaction 
(P < 0.0001) between orientation category and polar angles. One-way ANOVA and 
multiple-comparison tests indicated that horizontal and vertical orientations were 
over-represented in horizontal- and vertical-meridians {P < 0.0005 and P < 0.005, 
respectively). 



between preferred orientations and polar angles into 4 cat- 
egories (Radial, Radial +45, Tangential, and Radial— 45, see 
Fig. 6B,C). A 3-way ANOVA revealed significant main effects 
(P< 0.0001; P<0.05; P<0.01) of orientation category, the 
range of polar angles and the eccentricity, as well as a signifi- 
cant interaction (P< 0.0001) between orientation category and 
the range of polar angles. Subsequent 2-way ANOVA and 
multiple-comparison tests for each of the 3 ranges of polar 
angles indicated that the sensitivities for horizontal orien- 
tations were higher than other orientations at the horizontal 
meridian (P<0.01) and vertical orientations were higher than 
other orientations at the vertical meridian (P<0.05), respect- 
ively. In the oblique meridians away from horizontal and ver- 
tical meridians, we found that the sensitivities for both 
horizontal and vertical orientations (Radial —45° and Radial 
+45°) were significantly higher than those for oblique orien- 
tations (Radial and Tangential; P< 0.0001). The sensitivity for 
Radial orientations at oblique meridians appeared to be 
slightly higher than that for Tangential orientations, but this 
difference was statistically insignificant CP =0.21). These be- 
havioral results, resembling those reported previously (cf. 
Davey and Zanker 1998, their Fig. 1), thus, are grossly in 
agreement with the physiological findings revealed in our 
fMRI experiment, indicating the existence of perceptual 
superiority for cardinal orientations throughout the visual 
field, and in particular, near the horizontal and vertical 
meridians. 



Discussion 

Using high-resolution fMRI with a continuous and periodic 
stimulation paradigm and a differential analysis method, we 
reliably revealed orientation-tuning properties of single 
voxels in the human VI. At the population level, there was an 
overall trend for the voxels retinotopically located near 
horizontal- and vertical-meridian representations to respect- 
ively prefer horizontal and vertical orientations, thus support- 
ing and refining both the oblique effect and radial bias. 

Studying Single Voxels' Orientation Selectivity Using the 
Continuous and Periodic Stimulation Paradigm 

The continuous and periodic stimulation paradigm offers the 
advantages of exhaustive coverage of a given stimulus space 
and effective separation of the stimulus-evoked response from 
an intrinsic noise (Engel et al. 1994; Kalatsky and Stryker 
2003; Wandell et al. 2007; Engel 2011). In our initial 
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Figure 6. Discrimination thresholds of line orientations measured from 5 of the subjects who participated in the fMRI experiment. (4) Schematic drawing depicting 
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and C indicate 1 SD. 
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explorations, however, we found it difficult to reliably reveal 
single voxels' orientation preferences by using the original 
paradigm. As time went by, the BOLD response to the rotating 
grating either became weaker or lost its phase consistence 
with the stimulus orientation change. Two physiological 
mechanisms may account for this gradual deterioration. First, 
it is well documented that adaptation occurs in neural 
response (Movshon and Lennie 1979; Hammond et al. 1985) 
during prolonged stimulation, which is observable a few 
seconds after the stimulation onset. Secondly, the neurovascu- 
lar coupling, through which the BOLD signal reflects under- 
lying neural activity, may depend on the length of 
stimulation. Although it was initially believed that there is a 
tight coupling during prolonged stimulation (Bandettini et al. 
1997), recent results (Mintun et al. 2002; Lin et al. 2009) indi- 
cate that the coupling may be time-dependent; the coupling 
in the latter period (5-6 min after the stimulation onset) may 
not be as tight as that at the beginning. In the current study, 
we found a simple and effective solution to these 2 problems. 
By repeatedly presenting a full range of stimulus orientations 
and sandwiching these presentations in between extended 
blank periods, our paradigm takes the advantage of covering 
a complete stimulus space over a relatively short period 
(several tens of seconds) and minimizing the problems of 
adaptation and possible loose neurovascular coupling. 

Furthermore, as exemplified in Figure 1A, in comparison to 
the large elevation in the BOLD response to the grating, the 
fraction of the response modulated by the change in stimulus 
orientation was much smaller. The large elevation likely re- 
sulted from the mere presence of the grating and the spread 
of neurovascular responses. We then introduced a pair of 
orthogonal conditions, in which the identical grating was pre- 
sented with 2 different phases, one starting its rotation from 
horizontal orientation and the other from vertical orientation. 
The 2 orthogonal rotating gratings, when compared directly, 
thus, created a stimulus condition where the 2 sets of orthog- 
onal orientations changed continuously. Using the differential 
analysis method that makes a subtraction between the 
responses to the pair of orthogonal conditions, we removed 
the large elevation that was common to the 2 conditions and 
extracted the smaller response component that was modu- 
lated by the continuous change in stimulus orientation. 
Although the differential analysis method has been previously 
used in animal fMRI studies comparing a pair of orientations 
(e.g. 0° vs. 90° or 45° vs. 135°) with blocked designs (Duong 
et al. 2001; Fukuda et al. 2006; Moon et al. 2007; Serences 
et al. 2009), the present study provided the first evidence of 
its usefulness, when combined with blank periods, in study- 
ing orientation-tuning properties in the human visual cortex. 

Oblique Effect and Radial Bias in Orientation 
Representations in Human VI 

Orientation columns, as described for monkeys (Obermayer 
and Blasdel 1993), are spatially heterogeneous, consisting of 
several elements, including "linear zones", where orientation 
preferences change linearly, roughly 180° every 0.5-1.0 mm, 
"singularities," or point-like discontinuities, where orientation 
preferences change continuously through ±180° along any 
closed path around them, "fractures," where orientation pre- 
ferences change rapidly, and "saddle points," where orien- 
tation preference remains roughly constant. Studies in 



monkey VI using optical imaging have revealed close geo- 
metric relationships between the organizations of orientation 
preference and ocular dominance columns (Bartfeld and Grin- 
vald 1992; Obermayer and Blasdel 1993). Because ocular 
dominance columns in humans are approximately twice as 
wide as those in monkeys (Horton et al. 1990; Adams et al. 
2007), it may be assumed that the size of orientation hyper- 
columns in humans doubles that estimated for monkeys 
(0.5-1.0 mm, Blasdel 1992; Obermayer and Blasdel 1993), 
that is, orientation preferences change roughly 180° for every 
1-2 mm in humans. Thus, it is expected that using fMRI with 
an in-plane resolution of 0.625 x 0.625 mm 2 , the voxels 
sampled around linear zones, which amount to ~48% of area 
in VI (Obermayer and Blasdel 1993, their Table 2), are likely 
to show orientation preferences. On the other hand, it would 
be difficult to reliably reveal orientation selectivity around 
other elements (occupying >50% of area in VI), such as singu- 
larities, where neurons contained within a single voxel can 
respond to all orientations. Across 6 subjects, we found that 
more than a third of voxels (35-3%) had significant ccs, and 
were therefore considered to be orientation selective. Consid- 
ering other factors, such as the point spread function of the 
EPI approach we used (gradient-echo recalled EPI), which is 
~1.7-2 mm (e.g. Shmuel et al. 2007), and intrinsic blurring in 
EPI images due to T2* blurring, 35-3% was a very reasonable 
proportion. Because the in-plane resolution used in our study 
was still 1-2 orders of magnitude lower than that typically 
used in optical imaging, we refrained from using these voxels 
to construct spatial maps of orientation selectivity, which 
would require a extensive interpolation between identified 
orientation-selective voxels (Yacoub et al. 2008). Instead, we 
focused our analyses on global distribution features of indi- 
vidual voxel's orientation preferences. 

At the population level, we found that there were more 
voxels showing preferences for cardinal orientations (Fig. 4(7). 
It has been proposed that the behavioral superiority in cardi- 
nal orientation sensitivity (the oblique effect) originates from 
the differences in population size of neurons in the visual 
cortex (Maffei and Campbell 1970). Indeed, this anisotropy 
has been demonstrated in different species of animals using 
single-unit recording (e.g. Mansfield 1974; Li et al. 2003) and 
optical imaging (e.g. Chapman and Bonhoeffer 1998; 
Coppola, White et al. 1998; Wang et al. 2003). Our observation 
is in line with these findings in animals and provides direct 
evidence for the over-representations of cardinal orientations 
in the human VI. Alternatively, it has also been suggested that 
the behavioral superiority is linked to the difference in the am- 
plitudes of neuronal responses tuned to cardinal and oblique 
orientations (Maffei and Campbell 1970; Mansfield and 
Ronner 1978; Bonds 1982; Furmanski and Engel 2000). In 
several fMRI studies, it is reported that the oblique effect is 
manifested by either stronger (Furmanski and Engel 2000) or 
reduced (Serences et al. 2009; Mannion et al. 2010; Swisher 
et al. 2010) responses to cardinal orientations. In our own 
data, we only observed slightly larger responses to cardinal 
orientations than oblique orientations (Fig. 3D), but this differ- 
ence was not significant. It should be emphasized that our 
results were not affected by the fact that we only analyzed 
slightly more than a third of the voxels that showed significant 
orientation preferences. In additional analyses performed 
using all voxels in VI, we confirmed that enhanced response 
to particular orientations was not observed in the time course 
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of BOLD response (see Supplementary Fig. S5) or the strength 
of BOLD response (see Supplementary Fig. S6). The difference 
between our results and previous findings, thus, may be 
caused by the differences in the stimulation paradigm, spatial 
resolution and analysis method. On the other hand, the pre- 
vious observation that BOLD responses to cardinal orien- 
tations are lower than those to other orientations is less 
straightforward as far as interpretations are concerned. The 
reduced response to horizontal orientation, for example, was 
suggested to be a normalization mechanism to account for the 
fact that more cortical neurons respond to cardinal orien- 
tations than oblique orientations (Mannion et al. 2010). Future 
studies are warranted to elucidate these differences. 

In addition to the oblique effect, a radial bias, which de- 
scribes a biased representation of preferred orientations at 
particular angular positions, has also been shown in human 
psychophysics (Rovamo et al. 1982; Davey and Zanker 1998; 
Westheimer 2003, 2005; Sasaki et al. 2006), animal neurophy- 
siology (Leventhal 1983; Leventhal et al. 1984; Schall et al. 
1986; Bauer and Dow 1989; Rodionova et al. 2004; Durand 
et al. 2007), and both human (Sasaki et al. 2006; Mannion 
et al. 2009, 2010) and monkey (Sasaki et al. 2006) fMRI 
studies. One proposed neural mechanism for the radial bias is 
that the size of sub-regions representing different orientations 
in orientation columns changes depending on the retinotopic 
polar angle (Schall et al. 1986). The results obtained in our 
study, as well as those observed in several animal single-unit 
recording studies (Leventhal 1983; Leventhal et al. 1984; 
Schall et al. 1986; Bauer and Dow 1989) and a recent human 
fMRI study (Freeman et al. 2011), are in line with this hypoth- 
esis. By showing that horizontal and vertical orientations are 
predominantly represented around horizontal- and vertical- 
meridian representations, we also provided the first evidence 
that the radial bias was particularly striking for the cardinal 
orientations, suggesting that these portions of the cortex, at 
least in VI, are heavily engaged in processing cardinal orien- 
tations that are prevalent in the environment (Coppola, 
Purves et al. 1998; Hansen and Essock 2004; Girshick et al. 
2011). Conversely, our results do not support an alternative 
mechanism, proposed based on the results obtained from 
fMRI studies (Sasaki et al. 2006; Mannion et al. 2009, 2010), 
that BOLD responses to the radial orientations are greater 
than those to the orientations orthogonal to them. Our results 
are consistent with the findings of Freeman et al. (2011), indi- 
cating that the decoding of stimulus orientations using multi- 
variate pattern analysis (MVPA) approaches (e.g. Haynes and 
Rees 2005; Kamitani and Tong 2005) may in part rely on the 
radial bias. Moreover, if a spatial heterogeneity in the radial 
bias indeed exists, as suggested by our results, it may be pre- 
dicted that the decoding performance using MVPA varies 
between test orientations and sub-regions of the visual field. 

Conclusion 

In the present high-resolution fMRI study, we have introduced 
a stimulation and analysis method for efficiently revealing 
individual voxel orientation preferences in human VI. The 
basic idea of the stimulation paradigm and analysis method 
can be easily extended to study other stimulus tuning proper- 
ties, such as motion-direction preferences in medial temporal 
area (Albright et al. 1984; Zimmermann et al. 2011) and face 
identity and profile preferences in the inferotemporal cortex 



(Tanaka 2003), as well as certain types of modulatory effects, 
such as the contextual contrast modulation in early visual cor- 
tices (Tajima et al. 2010). 
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